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PREFACE 

This  work  was  undertaken  to  review  and  update  previous  evalua- 
tions of  thermal  imaging  systems'  performance  in  the  two  main  infra- 
red atmospheric  windows:  3-5  pm  and  8-12  pm.  Many  of  the  previous 
evaluations  were  based  on  outdated  or  erroneous  models,  especially 
those  for  atmospheric  transmission.  The  present  analysis  employed 
the  latest  available  models,  and  addressed  the  Air  Force  scenarios 
of  viewing  tactical  ground  target  from  airborne  systems. 
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SECTION  I 
SUMMARY 


A review  of  previous  analyses  for  the  relative  performance  of 
short  and  long  wavelength  thermal  imaging  systems  (i.  e. , nominal 
3-5  versus  8-12  (tm),  indicates  some  serious  shortcomings  in  the 
ancillary  models  employed.  All  analyses  consider  horizontal  sea  level 
atmospheric  paths,  neglect  or  approximate  path  radiance,  and  assume 
only  blackbody  functions  for  target  and  terrain  radiances.  The  present 
study  treats  atmospheric  slant  paths  from  various  aircraft  altitudes 
and  ranges  with  the  latest  Lowtran  transmission  model  (Lowtran  in), 
and  a sophisticated  path  radiance  model  developed  by  Anding  (Radst), 
for  various  seasonal,  latitudinal,  and  aerosol  models.  Actual  radi- 
ance measurements  of  typical  targets  and  terrains  are  employed,  in 
addition  to  theoretically  predicted  radiances  for  various  target  and 
terrain  temperatures. 

The  results  of  this  parametric  analysis  indicate  that  for  the 
nominal  target  considered  here  (i.  e. , tank,  rear  aspect),  the  short 
wavelength  system  is  superior  to  the  long  wavelength  system  at  all 
ranges  for  the  Tropical  and  Midlatitude  Summer  models  at  most  of 
the  typical  operating  conditions.  This  superiority  increases  with 
range,  until  at  twenty  kilometers,  there  is  about  an  order  of  magni- 
tude difference  between  the  performance  of  the  two  systems  for 
Summer  conditions,  and  about  two  orders  of  magnitude  for  Tropical 
conditions.  It  is  recognized  that  these  results  are  in  apparent  con- 
flict with  previous  analyses.  However,  comparisons  with  these  an- 
alyses do  indicate  general  agreement,  when  differences  in  assumed 
conditions  are  accounted  for.  For  Winter  conditions,  the  perform- 
ance of  the  two  systems  is  comparable  within  typically  30%  for 
ranges  from  zero  to  twenty  kilometers;  the  short  wavelength  system 
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is  slightly  better  at  shorter  ranges  and  slightly  worse  at  longer  ranges. 
For  cooler  targets  (or  at  other  aspects),  the  performance  of  the  short 
wavelength  systems  becomes  inferior  to  the  long  wavelength  systems 
at  short  ranges  for  the  Summer  and  Tropical  conditions,  and  at  all 
ranges  for  the  Winter  condition. 

Aerosols  are  predicted  to  have  little  effect  on  relative  perform- 
ance, but  this  is  thought  to  be  an  erroneous  result  due  to  a simplifi- 
cation in  the  Lowtran  model.  The  relative  performance  of  the  short 
wavelength  system  improves  with  increasing  target  temperature,  de- 
creasing altitude,  and  decreasing  terrain  temperature,  although  the 
latter  effect  is  not  particularly  significant.  No  large  advantages  are 
evident  for  other  spectral  bandpasses  within  these  two  window  regions, 
at  ranges  less  than  twenty  kilometers.  Comparisons  with  other  analyses 
of  f'ns  type  show  general  agreement,  considering  the  differences  in  con- 
ditions assumed.  In  order  to  evaluate  the  present  results  in  a laboratory 
environment,  an  infrared  transmissometer  and  peripheral  equipment  was 
delivered  in  August  1976,  for  use  in  the  AFAL  Tower  Facility  as  a cali- 
bration tool. 


SECTION  II 
INTRODUCTION 


Thermal  imaging  systems  (sometimes  referred  to  as  FLIRs, 
Forward  Looking  Mrared)  have  proven  their  value  for  a variety  of 
Air  Force  missions.  Most  of  these  systems  operate  in  the  8 to  12  pm 
atmospheric  window,  since  early  design  studies  indicated  this  window 
provided  superior  performance  over  the  3-5  pm  window  for  most  op- 
erating conditions.  More  recent  analyses  indicated  that  short  wave- 
length systems  may  be  superior  to  long  wavelength  systems,  for  long, 
humid  slant  paths.  There  are  also  some  operational  considerations 
which  may  favor  the  short  wavelength  system  (e.  g. , better  window 
thermal  properties,  better  resolution,  etc. ).  It  was  therefore  de- 
cided to  conduct  this  study  to  determine,  in  a parametric  fashion, 
the  situations  where  each  system  is  preferred.  Another  reason  for 
a new  analysis  is  that  many  of  the  previous  studies  were  based  on 
unsophisticated  or  outdated  models. 

The  scope  of  the  present  study  is  limited  to  an  aircraft  based 
system,  viewing  a typical  ground  target  at  ranges  less  than  twenty 
kilometers  from  altitudes  between  one-half  and  one  kilometer.  The 
latest  Lowtran  transmission  model  (i.  e. , Lowtran  HI)  was  to  be  used 
to  consider  various  atmospheric  models.  The  range  of  parameters 
employed  in  this  study  is  outlined  in  Table  1.  In  general,  the  an- 
alysis concentrated  on  using  updated  atmospheric,  target  and  back- 
ground models,  while  using  an  existing  sensor  performance  model. 

It  should  be  emphasized  that  this  study  is  limited  to  analyzing  the  per- 
formance of  a short  wavelength  infrared  (i.  e. , SWIR,  defined  here 
nominally  as  3.4  - 4.8  pm)  thermal  imaging  system  relative  to  a long 
wavelength  infrared  (i.  e. , LWIR,  defined  here  nominally  as  8. 1 - 
12.2  pm)  system.  Whether  either  system  will  operate  adequately  for 
the  conditions  considered  here  is  not  within  the  scope  of  this  study. 
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Table  1. 

Range  of  Study  Parameters 


[ 


Platform 

Altitude 

Range 

Sensor 

Spectral  Region  s* 
Atmospheres 

Lowtran  in  Models 

Aerosol  Model 
Continuum 

Target  (Typical  Tank) 

Aspect 

Temperature 

Background 

Terrain  Radiance 

Path  Radiance 


0.5,  0.75*,  1.0  km 
0-20  km 

3.4  - 4.8  and  8.1  - 12.2  pm 


Midlatitude  Winter  and  Summer^  and 
Tropical 

5,  8.5*,  23  km  Visibility 
Burch  4 pm  Water  Continuum2* 


Front,  Side,  Rear* 

Measured  and  250,  275,  325,  350  K 


Day  and  Night*  Measured  and 
250,  270,  290,  310  K 

Same  Conditions  as  "Atmospheres,  " Above  ' 


* Denotes  nominal  condition  which  were  generally  employed  when  the 
other  parameters  were  varied. 

* Other  narrower  bands  within  these  nominal  bands  were  also  in- 
vestigated. 

2 Other  continua  at  4 and  10  pm  were  also  investigated. 


The  organization  of  this  report  is  as  follows.  Major  conclusions 
and  recommendations  are  summarized  in  Section  L The  purpose  and 
scope  of  the  study,  and  the  organization  of  the  report  are  given  in  the 
present  section  (i.  e.,  Section  II).  Section  III  surveys  some  of  the  more 
important  existing  literature  in  this  field.  The  theoretical  development 
and  main  characteristics  of  the  various  models  employed  in  the  present 
analysis  (i.  e. , the  sensor  model,  the  atmospheric  transmission  and 


f 
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radiance  models,  and  the  target  and  terrain  radiance  models)  are  given 
in  Section  IV.  Calculated  relative  performance  results,  for  thermal 
imaging  systems  operating  in  the  short  and  long  wavelength  atmospheric 
windows,  are  presented  in  Section  V.  Finally,  Section  VI  indicates  the 
main  conclusions  of  the  analysis,  and  makes  recommendations  for  future 
investigations  in  this  field.  Ancillary  information  is  also  presented  in 
three  appendices  (Appendix  A,  B and  C). 
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SECTION  III 
LITERATURE  SURVEY 


Biberman  considered  the  problem  of  detecting  a target  against 
uniform  and  cluttered  backgrounds  [1].  He  proposed  the  following 
criteria: 

2 

J/R  ^ u>L  uniform  background 

2 (3-1) 

J/R  ^ ^u^raax  cluttered  background 


where  J is  the  target  radiant  intensity,  R is  the  range,  to  is  the  field 
of  view,  and  L is  the  background  radiance.  For  the  cluttered  case, 
the  criterion  is  necessary  but  perhaps  not  sufficient,  and  the  maximum 
value  to  background  radiance  in  the  field  of  view  must  be  used.  This 
discussion  did  not  consider  such  parameters  as  optical  gain,  cell  noise, 
cell  sensitivity  or  system  noise,  and  did  not  give  any  quantitative  results. 


Mundie  considered  the  relative  merits  of  detectors  which  utilize  the 
3-5  and  8-13  /im  atmospheric  windows  [2].  He  suggested  the  following 
performance  criteria: 


Signal-to-Noise  Ratio: 


Clutter-to-Noise  Ratio: 


Contrast: 


C 


S . -B  . 
min  mm 

L 

B - B . 

max  mm 

L 

S . -B 
= nun  max 

B -B  . 

max  mm 


(3-2) 


Detectability: 


5 . -E 

™ = mm  max 
N 


where  B . and  B are  the  averages  of  the  detector  output  sig- 
mm  max 

nals  generated  by  viewing  members  of  the  hottest  and  coolest  background 
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families,  respectively;  and  Smin  and  Smax  are  the  signals  when  a 
small  target  is  in  the  field  of  view  against  the  respective  background 
families.  Using  these  expressions,  Mundie  evaluated  the  relative 
performance  of  short -to-long  wavelength  FLIRs  assuming  blackbody 
radiances  for  the  targets  and  backgrounds,  and  a "temperate”  at- 
mospheric transmission  model.  The  transmission  was  representative 
of  50  atm  • cm  of  CC>2  and  2 pr  • cm  of  HgO,  based  on  measurements 
by  Yates  and  Taylor  [3],  and  Wyatt,  et.  al.  [4,  5].  He  concluded  that 
the  performance  of  the  short  wavelength  scanner  was  inferior  to  the  long 
wavelength  scanner  by  factors  of  2.5  and  8.4  at  ranges  of  zero  and  1.6  km, 
respectively.  The  calculated  performance  of  the  two  systems  is  equiv- 
alent at  target  temperatures  of  440  and  530  K,  at  ranges  of  zero  and 
1.6  km,  respectively. 


Rosell,  Sendall,  et.  al.,  considered  the  image  contrast  and  modu- 
lation contrast  for  uniform  and  variegated  oackgrounds,  respectively 
[6-14]: 


-L  . 

max  mm 


max 


L -L  . 

max  min 

L + L . 

max  mm 


uniform  background 


variegated  background 


(3-3) 


where  L is  the  radiance  of  the  target  or  background.  They  developed 


complex  expressions  for  system  performance  in  uniform  and  cluttered 
backgrounds  considering  such  parameters  as  finite  apertures  of  the 
sensors  and  other  elements  of  the  system  (i.  e. , modulation  transfer 
functions),  the  image  size  and  shape,  the  integration  time  of  the  eye, 
picture  aspect  ratio,  etc.  However,  they  apparently  did  not  carry 
out  calculations  on  relative  performance  of  various  systems. 
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Schnitzler  considered  composite  FLIR-visual  systems'  perform- 
ance in  detecting  targets  of  various  sizes  [15, 16].  He  developed  ex- 
pressions for  the  "noise  required-input  temperature  difference" 

At.  to  express  the  performance  of  thermal  imaging  systems  for  spot 
detection.  Using  values  of  the  various  psychophysical  and  system  param- 
eters based  on  test  data,  the  following  expressions  resulted: 


At.  =7.19x10 
in 


AT.  = 6.83  x 10 
m 


-MS)] 

-•MS)] 


8.13  - 12.2  *im 


3.33  - 4.17  pm 


(3-4) 


where  is  the  background  radiance  and  z is  a function  of  the  angular 
width  of  the  target,  a. 

Using  Lowtran  II  with  the  standard  Tropical  atmospheric  model  at 
5 and  23  km  visibility  [17],  ATin  was  evaluated  as  a function  of  both 
a and  range.  According  to  these  calculations,  the  long  and  short  wave- 
length FLIR  systems  give  equivalent  performances  between  about  3.5 
and  3.9  km  for  both  visibilities  and  all  target  sizes;  the  short  wave- 
length system  is  superior  at  larger  ranges. 

Barhydt,  et.  al. , developed  a "radiation  function"  to  compare 
long  and  short  wavelength  thermal  imaging  infrared  sensors  [18, 19]: 


J Tf(X)Ta(X)(lAVd> 

1 0 1 a \6T&V 


2Vr5hc" 


“ t>L* 

jFiW-bt* 


where  t,  and  r are  the  transmissions  of  the  filter  and  the  atmosphere, 

f 3. 

respectively,  Jt  is  the  radiant  intensity  of  the  target,  and  L*  is  the 
radiance  of  the  background  in  quantum  units  (e.  g. , photons/sec/sr/m  ). 
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Other  spectral  dependent  parameters  are  also  considered  to  evaluate 
the  relative  performance: 

RAT  = ^d^q  RAT 

where  K is  unity  for  a photoconductive  detector  and  for  a photo- 
voltaic detector,  is  the  detector  quantum  efficiency,  and  is  the 
detector  size  parameter.  Four  measurements  of  atmospheric  trans- 
mission by  Yates  and  Taylor  [3],  at  different  ranges  and  humidities, 
were  employed  to  evaluate  this  function;  300  K blackbodios  were  assumed 
for  the  target  (J^)  and  the  background  (L* ).  For  two  different  system 
aperture-resolution  combinations,  the  performance  of  the  short  wave- 
length system  was  calculated  to  be  1.7  to  2.7  times  better  than  the  long 
wavelength  system  at  a range  of  18.7  km  and  a relative  humidity  of 
82%  (see  Table  2).  At  the  other  three  range-humidity  combinations, 
the  short  wavelength  system  was  calculated  to  be  inferior  to  long  wave- 
length system. 


Table  2. 

Relative  Performance  from  Barhydt  [19]; 
3.4-4.8  vs  8.1-12.2  pm  Systems 

Range 

(km) 

Relative 

Humidity 

<%) 

Aperture - 
Resolution 
(in-mrad) 

Relative 

Performance 

5.5 

51 

1.  25 

0.  30 

5.  5 

73 

1.25 

0.  33 

16.3 

53 

1.  25 

0.  37 

16.3 

82 

1.  25 

1.  74 

5.5 

51 

0.8 

0.  46 

5.5 

73 

0.8 

0.51 

16.3 

53 

0.8 

0.  57 

16.3 

82 

0.8 

2.74 

Milton,  et.  al.  [20],  used  the  performance  model  of  Barhydt  [18] 
in  a similar  analysis,  except  they  used  the  Lowtran  II  atmospheric 
transmission  model  [17]  to  evaluate  r , instead  of  using  the  measure- 
ments  of  Yates  and  Taylor.  Calculations  were  presented  as  a function 
of  range,  for  three  visibilities  (i.  e. , aerosol  models)  and  four  humidi- 
ties.  At  a humidity  of  18  gm/m  (near  the  Tropical  model  Schnitzler 
used),  the  long  and  short  wavelength  systems  were  predicted  to  give 
equivalent  performances  at  ranges  from  about  5 to  7 km. 

Moser  [21]  referred  to  a mathematical  model  of  airborne  FLIR 

operational  performance,  but  did  not  describe  it.  He  did  reference 

several  reports  [22-24],  but  due  to  time  limitations  of  the  current 

contract,  these  reports  were  not  reviewed.  Shumaker  and  Keller  [25] 

presented  a model  of  FLIR  performance,  but  did  not  show  any  com- 

* 

parisons. 

In  summary,  the  four  authors  who  presented  relative  performance 
were  Mundie,  Schnitzler,  Barhydt  and  Milton  (their  results  are  com- 
pared to  the  present  results  in  Section  V).  All  used  simple  blackbody 
radiances  for  the  target  and  background.  For  atmospheric  transmis- 
sion values,  Mundie  and  Barhydt  used  measurements  by  Yates  and 
Taylor,  while  Schnitzler  and  Milton  employed  the  Lowtran  II  code. 

All  authors  considered  only  horizontal  sea-level  paths,  and  apparently 
neglected  the  path  radiance  contribution  to  background  noise  (see  Sec- 
tion IV,  #3). 


SECTION  IV 

THEORETICAL  DEVELOPMENT 


i 


In  this  section  formulations  are  given  for  various  equations  and 
procedures  used  in  the  calculation  of  the  relative  performance  of  ther- 
mal imaging  systems.  Section  IV,  #1  gives  the  development  of  the 
performance  model,  Section  IV,  #2  discusses  the  atmospheric  trans- 
mission considerations,  Section  IV,  #3  considers  background  radiance, 
and  Section  IV,  #4  presents  target  signature  calculations. 


1.  Performance  Model 

A model  developed  by  Barhydt  [18]  was  used  in  this  study  to  eval- 
uate the  relative  performance  of  thermal  imaging  systems: 


b2J, 


t 


R 


AT 


2 J 2hc 


rf  (X)  Ta(X)  XdX 
t>L£  ~|l/2 


(4-D 


For  relative  performance,  only  spectrally  varying  factors  are  impor- 
tant, so  the  constant  in  front  of  the  equation  can  be  neglected.  The 
filter  function  (r^)  is  assumed  to  be  perfect,  i.  e. , unity  inside  the 
spectral  band  and  zero  outside.  The  X in  the  numerator  is  simply  to 
correct  for  not  using  quantum  units  in  the  numerator  for  J.  If  photons 
are  used  here  (i.  e. , J*),  as  in  the  numerator,  X can  be  removed. 
Thus  we  have: 


R 


J 

^X 


b2jt 

Ta^bTbX  dX 


at  f bL*  ni/2 


(4-2) 
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Three  other  spectrally  dependent  terms  from  the  full  signal-to-noise 
expression  must  also  be  considered: 


rat  = kV\rat  (4-3) 

where  K is  unity  for  photoconductive  detectors,  but  equal  to  J~2  for 
photovoltaic  detectors;  V is  the  detector  quantum  efficiency;  and 

4 

is  the  normalized  detector  size  parameter.  Barhydt  [19]  recom- 
mends values  of  0.5  for  v for  detectors  in  both  spectral  regions* 

(e.  g. , indium  antimonide  in  the  3-5  (im  region,  and  mercury-doped 
germanium  in  the  8-12  pm  region).  To  compensate  for  the  inherently 
lower  value  of  MTF  in  the  8-14  Mm  region  (due  to  the  increase  in  op- 
tical diffraction,  smaller  detectors  can  be  employed.  To  obtain  the  same 
MTF,  T7d  for  long  wavelength  systems  is  reduced  from  unity  in  the  short 
wavelengths  to  between  0.5  and  0.8  (for  aperture  diameter  - nominal 
angular  resolution  products:  Dp  = 0.8  - 1.25  in-mrad,  respectively 
[18]).  Thus,  the  short  to  long  wavelength  ratio  of  these  three  param- 
eters (i.  e. , the  relative  detector  performance,  £)  varies  between: 


t*  (K^>”  . ,.77  - 2.83 
(k^)lwir 


(4-4) 


In  the  present  study,  the  larger  value  (i.  e. , 2.83  corresponding  to  more 
advanced  systems)  was  used.  **  Finally,  the  relative  performance  is 
expressed  as: 


P = 4 


(ratWir 

^at^lwir 


(4-5) 


* In  his  original  paper  [18],  Barhydt  suggested  tj  = 0.25  for 
mercury-doped  germanium  detectors.  ^ 

**  It  should  be  noted  that  neither  Mundie  [2]  or  Milton  [20]  included 
this  factor. 
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2.  Atmospheric  Transmission 

Previous  studies  of  this  nature  (see  Section  III)  have  used  either 
measurements  by  Yates  and  Taylor  [3]  or  the  Lowtran  II  code  [17] 
for  evaluating  the  spectral  atmospheric  transmission.  Although  mea- 
surements (such  as  those  of  Yates  and  Taylor)  are  probably  more 
realistic  than  calculations,  they  are  limited  in  the  range  of  conditions 
which  can  be  considered.  Lowtran  is  a one  parameter  band  model 
for  calculating  atmospheric  transmission  for  arbitrary  slant  paths 
based  on  techniques  suggested  by  Altshuler  [26].  It  also  considers 
aerosol  extinction  (i.  e. , scattering  and  absorption),  and  certain  con- 
tinuum molecular  absorption. 

a.  Molecular  Line  Absorption 

The  Lowtran  II  model,  used  in  some  previous  studies  of  this  type, 
has  been  superceded  by  Lowtran  III  which  updates  some  of  the  band 
model  parameters  [27].  Lowtran  considers  the  uniformly  mixed  at- 
mospheric gases:  CC>2,  N20,  CH4,  CO,  and  02;  as  well  as  non- 
uniform  water  vapor  and  ozone.  Its  main  assumption  is  that  the  aver- 
age transmittance  over  a twenty  wavenumber  spectral  interval  can  be 
represented  by  a single  parameter  of  the  form: 

f = f(C„w)  (4_6) 

where  is  the  wavelength  dependent  absorption  coefficient  and  w is 
the  equivalent  absorber  amount.  Lowtran  is  most  inaccurate  in  the 
window  regions  between  strong  bands  (i.  e. , the  weak  line  approxima- 
tion region).  A warning  is  given  [27]:  "The  one-dimensional  predic- 
tion scheme  presented. . . is  not  accurate  for  such  conditions. . . ". 

These  of  course  are  iust  the  regions  of  most  interest  in  the  present 
study. 
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A sample  line-by-line  calculation  (i.  e. , Hitran)  using  the  AFCRL 
line  parameter  compilation  [28],  was  performed  for  comparison  with 
the  Lowtran  in  results.  Because  of  the  expense  of  such  line-by-line 
calculations,  only  two  narrow  regions  in  the  short  wavelength  window 
(i.  e. , 3.4 -4.1  and  4.55-4.  8 pm)  were  considered.  To  be  consistent 
with  the  Lowtran  results,  the  infinite  resolution  line-by-line  results 
were  degraded  to  20  wavenumber  resolution  (using  a square  filter 
function).  Figures  1-4  compare  the  spectral  transmission  calculated 
by  the  two  techniques  for  the  Midlatitude  Summer,  sea  level  condition 
at  several  ranges.  It  is  immediately  apparent  that  although  the  re- 
sults were  presumably  degraded  to  the  same  resolution,  the  line-by- 
line calculations  show  much  more  spectral  detail  than  those  of  Low- 
tran. Note  also  that  the  degree  of  agreement  at  various  wavelengths 
is  different  at  different  ranges  (some  improving  and  some  becoming 
worse  with  range).  This  is  characteristic  of  a single  parameter  band 
model,  and  can  be  improved  using  a two-parameter  model. 

A comparison  of  the  transmission  integrated  over  the  spectral 
band  is  shown  in  Figure  5 as  a function  for  range.  It  is  seen  that 
for  the  3.4  -4.1  pm  spectral  region,  Lowtran  gives  low  transmission 
relative  to  the  line-by-line  calculations  at  short  ranges,  and  high 
transmission  at  longer  ranges.  The  maximum  error  for  this  band 
for  ranges  less  than  40  km  is  about  ± 6%.  For  the  other  bandpass 
(i.  e.  , 4.55  -4.8  pm),  the  error  is  about  the  same  at  short  ranges, 
but  increases  with  range  beyond  about  20  km.  At  40  km  range,  the 
Lowtran  results  are  almost  a factor  of  two  below  the  line-by-line 
results.  Water  continuum  and  aerosol  extinction  are  not  included  in 
these  results  (nitrogen  continuum  is,  since  it  cannot  easily  be  sup- 
pressed from  the  Lowtran  results).  However  since  these  components 
(i.  e.  , water  continuum  and  aerosol  extinction)  are  slowly  varying 
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Figure  1 . Comparison  of  Lowtran  and  Hitran  Calculated  Atmospheric  Spectral  Transmission 
Including  Nitrogen  Continuum,  Midlatitude  Summer  Sea  Level  Model,  1 km  Range 


e 2.  Comparison  of  Lowtran  and  Hitran  Calculated  Atmospheric  Spectral  Transmission 
Including. Nitrogen  Continuum,  Midlatitude  Summer  Sea  Level  Model,  10  km  Range 


Figure  3.  Comparison  of  Lowtran  and  Hitran  Calculated  Atmospheric  Spectral  Transmission 
Including  Nitrogen  Continuum,  Midlatitude  Summer  Sea  Level  Model,  20  km  Range 
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Figure  4.  Comparison  of  Lowtran  and  Hitran  Calculated  Atmospheric  Spectral  Transmission 
Including  Nitrogen  Continuum,  Midlatitude  Summer  Sea  Level  Model,  40  km  Range 


functions  of  wavelengths  (see  Sections  IV,  lb  and  IV,  lc),  they  are,  to 
a good  approximation,  multiplicative  factors.  Thus,  the  ratio  shown 
in  Figure  5 should  be  approximately  unchanged  with  the  inclusion  of 
those  attenuations. 

b.  Molecular  Continuum  Absorption 

In  the  window  regions,  there  is  another  important  molecular  ab- 
sorption mechanism:  continuum  absorption.  In  the  short  wavelength 
region  carbon  dioxide,  nitrogen  and  water  absorb  in  a continuous 
fashion.  Carbon  dioxide  absorption  is  small  compared  to  other  mech- 
anisms [29],  and  is  not  included  in  Lowtran  or  in  the  present  study. 

Nitrogen  modeling  in  Lowtran  is  based  on  measurements  of  Reddy 
and  Cho  [30],  and  Shapiro  and  Gush  [31].  Other  measurements  by 
Burch  [29]  are  more  widely  used,  having  been  employed  by  Spencer 
[32],  Long  [33]  and  Meredith  [34]  in  their  modeling.  Burch's  mea- 
surements are  in  general  agreement  with  the  Lowtran  values  (Low- 
tran is  only  about  20%  less  absorbing  at  4 /im).  Thus  for  convenience, 
the  existing  Lowtran  nitrogen  continuum  was  used  in  the  present  study 
(see  Figure  6). 


Figure  6.  Lowtran  Nitrogen  Continuum  Transmission 
Model;  5 km  Range,  Midlatitude  Summer 
Atmospheric  Model  at  Sea  Level 
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The  Lowtran  III  model  was  also  supposed  to  include  the  4 *xm 
water  continuum  [27]:  "Absorption  coefficients  for  the  water  con- 
tinuum near  10  pm  and  4 pm  in  Lowtran  III  are  based  on  measure- 
ments. . . However,  in  examining  the  output  and  the  code,  it  is 
apparent  that  this  continuum  is  still  being  neglected. 

There  is  still  some  debate  about  the  proper  level  of  the  contin- 
uum near  4 pm.  The  generally  accepted  level  is  based  on  extrapola- 
tions of  data  taken  by  Burch  at  high  temperature,  high  pressure,  and 
low  spectral  resolution  [27].  Because  the  continuum  absorption  is 
small  and  has  little  spectral  structure,  measurements  are  difficult 
even  with  laboratory  path  lengths  of  a kilometer.  Recently  however, 
two  sets  of  White  cell  measurements  of  total  water  vapor  absorption 
coefficients  have  been  reported,  which  correspond  to  Midtatitude 
Summer  conditions  [35,  36].  In  general,  it  seems  that  currently  used 
extrapolations  from  Burch  data  are  somewhat  higher  in  transmission 
than  reality.  However,  they  are  felt  to  be  valid  to  within  a factor  of 
two  [37].  This  continuum  was  incorporated  into  the  Lowtran  results 
for  the  present  study  using  a special  purpose  computer  code  "Tran" 
(see  Appendix  A). 

In  the  long  wavelength  window,  the  only  important  continuum  ab- 
sorption is  due  to  water.  There  is  considerable  controversy  regarding 
this  continuum.  It  is  generally  evaluated  from: 

k “ ww  (Cspw  + Cfpf)  <4‘7) 

where  k is  the  continuum  absorption  coefficient,  ww  is  the  water 
vapor  content,  pw  and  p{  are  the  partial  pressure  of  water  vapor 
and  foreign  gases,  and  Cg  and  Cf  are  the  self  and  foreign  broadening 
absorption  coefficients.  Lowtran  uses  measurements  of  Burch,  et.  al. 
[38,  39],  McCoy  and  Rensch  [40],  and  Bignell  [41],  for  the  self  broad- 
ening absorption  coefficient,  Cg. 
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Two  significant  approximations  are  then  invoked: 


* 


I 

** 
• . 


1.  The  self  broadening  absorption  coefficient  (C  ) does  not 

b 

vary  with  temperature. 

2.  The  foreign  broadening  absorption  coefficient  (Cf)  is 

proportional  to  the  self  coefficient  (C  ). 

s 

The  first  approximation  was  probably  used  for  lack  of  accurate  know- 
ledge of  the  temperature  effect.  The  second  was  based  on  a single 
measurement  by  McCoy,  et.  al.  [40],  with  large  experimental  un- 
certainty [37].  The  value  used  for  the  ratio  C./C  is  0.005. 

I s 

Roberts,  et.  al.  [42]  have  proposed  a modified  expression  which 
lifts  the  first  approximation,  but  not  the  second.  It  does  however, 
alter  the  value  of  Cj/Cg  to  0.0008,  based  on  more  recent  measure- 
ments by  Long  and  Burch  [42].  The  level  of  C was  also  modified 

s 

slightly  based  on  more  recent  measurements  by  Eureh,  The  tempera- 
ture dependence  of  C was  taken  to  be: 

s 

Cs=c^[to<¥'255>]  (4-8) 

where  C°  is  the  value  of  C3  at  296  K,  and  T is  the  temperature  in 
Kelvin.  Roberts  suggests  a value  of  1800  K for  T . 

It  should  be  noted  that  due  to  the  second  approximation,  has 

the  same  temperature  dependency  in  this  model  as  C . There  is  no 

s 

reason  to  expect  the  same  dependency,  since  two  different  mechanisms 
are  involved.  Indeed,  there  is  evidence  that  the  two  parameters  should 
have  opposite  temperature  effects.  A model  by  Kunde  [43],  treats  C 

s 

and  Cj  separately,  and  assigns  them  different  (and  opposite)  tempera- 
ture dependence. 

Sample  results  for  these  three  models  are  shown  in  Figure  7. 

It  is  seen  that  Roberts'  model  varies  with  temperatui  n much  more 
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7.  Comparison  of  Various  8-12  pm  Water  Continuum  Models. 
Water  content  (14  gm/m^)  and  total  pressure  (1  atm) 
correspond  to  Midlatitude  Summer  Model. 
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Transmission  (%)  (10  km  Path) 


than  the  others  because  both  C and  C.  are  assumed  to  have  the  same 

O I 

functional  behavior.  Although  C (and  C,)  in  the  current  Lowtran 

S I 

model  are  constant  with  temperature,  k varies  with  temperature  be- 
cause the  water  vapor  partial  pressure  (p  ) is  proportional  to  tern- 
perature  for  a given  water  content  w . It  is  felt  that  the  current 
Lowtran  model  for  10  pm  water  vapor  continuum  is  in  error,  but 
that  Roberts'  model  probably  over-corrects  for  the  temperature  effect. 
Kunde’s  model  is  probably  the  best  choice,  but  since  there  is  no  strong 
evidence  supporting  it,  the  existing  Lowtran  m model  was  employed 
in  this  study.  It  will  be  shown  that  the  various  continuum  models  are 
not  too  far  apart  for  the  conditions  considered  here. 


c.  Aerosol  Extinction 

Aerosol  extinction  is  the  result  of  both  scattering  and  absorption, 

and  depends  on  the  size  distribution  of  the  aerosol  particles  ^ and 

dr 

on  their  optical  properties  (i.  e. , complex  index  of  refraction,  n): 


kext~*f  Qflvf^»n»r)  ^Sr^^dr 


ext 


dr 


(4-9) 


where  Qexj.  is  the  aerosol  extinction  efficiency  function,  X is  wave- 
length and  r is  the  radius  of  the  particle.  The  particle  size  distribu- 
tion is  probably  the  most  important,  and  most  variable  factor  affecting 
aerosol  extinction.  Figure  8 shows  the  variety  of  size  distributions 
which  have  been  used  by  various  investigators  in  attempting  to  des- 
cribe various  real  world  aerosol  size  distributions  with  meteorological 
and  geographical  parameters.  Variables  which  are  often  included  in 
such  models  are  relative  humidity,  wind  speed,  and  the  influence  of 
oceanic  or  urban  environments. 


The  complex  index  of  refraction  of  the  aerosol  affects  the  aerosol 
extinction  through  the  efficiency  function,  Q . For  typical  aerosol 
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materials  (e.  g.  , dust,  soot,  sea  salts  and  water  solubles),  the  real 

part  of  the  index  has  values  in  a relatively  narrow  range  from  about 

1.2  to  1.8  for  the  4 and  10  pm  regions  [44].  The  imaginary  part  of 

the  index  for  these  materials  and  spectral  regions,  ranges  from 
-3 

about  10  to  almost  unity,  but  its  effect  on  the  extinction  coefficient 
is  smaller  than  that  of  the  real  index.  For  a typical  measured  size 
distribution,  calculations  indicate  that  the  aerosol  extinction  coeffi- 
cient varies  by  about  a factor  of  two  for  these  ranges  of  optical  prop- 
erties [45]. 

Several  models  are  available  for  predicting  aerosol  extinction. 
Figure  9 indicates  the  results  of  several  of  these  models  as  a func- 
tion of  relative  humidity  for  the  4 pm  region.  The  curves  in  this  fig- 
ure are  for  the  models  of  Gambling  [46],  Hodges  [47],  and  Barnhardt 
and  Streete  [48]  in  which  humidity  is  a parameter.  The  various  sym- 
bols in  the  figure  depict  the  models  of  Katz  [49],  AFCRL  [50-54],  and 
Hughes  [55].  The  numbers  in  the  parentheses  indicate  the  generic 
type  of  aerosol  (e.  g.  , (1)  - all  maritime  type  aerosol,  (0)  - all  con- 
tinental). Figure  10  shows  recent  NRL  measurements  of  aerosol 
extinction  near  4 pm  [56]  superimposed  on  the  models  given  in  Fig- 
ure 9.  Also  shown  in  this  figure  are  British  measurements  [57,  58], 
which  have  not  been  corrected  for  molecular  absorption.  These 
models  and  measurements  show  the  same  general  trend  with  relative 
humidity;  however,  the  variability  is  large  at  a given  humidity.  Sim- 
ilar results  are  shown  for  the  10  pm  region  in  Figure  11. 

Lowtran  III  uses  the  AFCRL  model  indicated  (clear  and  hazy, 
designated  in  Lowtran  by  "visibilities"  of  23  and  5 km,  respectively) 
in  Figures  9-11.  It  can  be  seen  that  the  extremes  in  extinction  re- 
sulting from  this  model,  almost  bracket  the  extremes  indicated  by 
the  other  models  and  measurements.  It  was  therefore  decided  to 
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Figure  9.  Comparison  of  Aerosol  Extinction  Models 


Figure  11.  Comparison  of  Aerosol  Extinction  Measurements 
and  Models  in  10  pm  Region 


employ  this  model  (with  one  intermediate  condition,  i.  e. , 8.5  km 
visibility)  in  the  present  calculations. 

A note  of  caution  regarding  the  Lowtran  aerosol  model  should 
be  interjected  here.  The  model  adjusts  only  the  total  aerosol  num- 
ber density  (see  Figure  8,  curves  f and  g)  to  arrive  at  different 
visibilities  [27],  but  does  not  adjust  the  relative  size  distribution 
(e.  g. , the  relative  number  of  larger  versus  smaller  particles),  nor 
the  optical  properties  of  the  aerosols.  Thus,  the  spectral  shape  of 
the  aerosol  extinction  is  similar  for  different  visibilities  (see  Fig- 
ure 12).  This  approximation  could  give  erroneous  indications  of  the 
effects  of  aerosols  on  the  relative  performance  of  thermal  imaging 
systems  operating  in  the  two  atmospheric  windows. 

d.  Sample  Results 

The  atmospheric  transmission  factor  required  to  evaluate  the 
performance  of  thermal  imaging  sensors  was  calculated  using  the 
Lowtran  in  code,  and  incorporating  the  4 pm  water  continuum. 
Sample  results  of  atmospheric  transmission  for  various  conditions 
are  shown  in  Figures  13-26.  Transmission  as  a function  of  range 
is  shown  in  Figures  13-18  for  nominal  atmospheric  models  (i.  e. , 
AFCRL  Midlatitude  Winter  and  Summer,  and  Tropical  atmospheres) 
with  8.5  km  visibility,  and  a nominal  platform  altitude  of  0.75  km. 
The  strong  feature  near  4.3  pm  is  due  to  carbon  dioxide  absorption, 
while  most  of  the  structure  between  3.5  and  3.9  pm  is  caused  by 
methane.  The  long  wavelength  curves  have  less  structure,  and  are 
generally  lower  than  the  short  wavelength  curves.  Most  of  the  ab- 
sorption in  this  region  is  caused  by  water  (both  line  and  continuum), 
with  some  ozone  and  carbon  dioxide  attenuation.  The  features  at  9.5 
and  9.7  pm  are  due  to  ozone,  but  most  of  the  other  features  are  from 
water.  The  cutoff  near  8 pm  is  due  to  methane  and  NgO.  Because 
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Figure  12.  Lowtran  III  Calculated  Aerosol  Spectral  Transmission  for  Two 
Visibilities  in  the  Long  Wavelength  and  Short  Wavelength 
Infrared  (LWIR  & SWIR)  Atmospheric  Windows 


ltitude 


Figure  18.  Calculated  LWIR  Spectral  Atmospheric  Transmission  as  a Function  of  Range; 
Tropical  Model  Atmosphere,  8.5  km  Visibility,  0.75  km  Altitude 


of  the  lower  water  content  in  the  Winter  model,  its  transmission  is 
much  better  than  the  Summer  model’s,  particularly  in  the  long  wave- 
length window  where  water  is  more  important. 

Figures  19  and  20  present  these  results  at  a particular  range  as 
a function  of  atmospheric  model,  mainly  to  indicate  the  effect  of  water 
content.  Long  wavelengths  are  considerably  more  affected  than  short 
wavelengths,  because  of  the  more  dominant  water  continuum  and  line 
absorption  in  that  region  (i.  e. , 8-12  pm).  It  is  interesting  that  some 
narrow  bands  inside  the  short  wavelength  region  (e.  g. , 3. 8 -4.1  pm) 
show  very  little  reduction  with  humidity.  These  may  be  good  candidate 
regions  for  systems  operating  in  very  humid  locales  (see  Section  V). 

The  effect  of  the  various  water  continua  models  (see  Section  IV,  #b), 
is  shown  in  Figures  21  and  22  for  a nominal  Tropical  condition.  For 
the  lower  water  content  of  the  Summer  or  Winter  conditions  the  varia- 
tions would  be  much  less.  It  is  seen  that  the  short  wavelength  contin- 
uum is  only  significant  between  about  3.5  and  4.2  pm.  For  the  long 
wavelength  region,  the  differences  between  the  various  continuum 
models  would  be  much  less  pronounced  at  the  Summer  condition  be- 
cause of  the  lower  temperature  (the  average  temperature  for  this  slant 
path  is  about  292  K,  while  the  cross-over  points  for  the  models  are 
around  283  K,  see  Figure  7),  as  well  as  the  reduced  water  content. 

I 

The  effect  of  various  aerosol  models  (i.  e. , visibility,  see  Sec- 
tion IV,  #c),  on  the  atmospheric  transmission  is  shown  in  Figures  23 
and  24  for  the  nominal  model  atmosphere  (i.  e. , Midlatitude  Summer) 
at  nominal  conditions.  It  is  seen  that  the  aerosol  extinction  is  some- 
what less  at  the  long  wavelengths  than  at  the  short.  The  extinction 
due  to  aerosols  is  assumed  independent  of  season  and  locale  (i.  e. , 
the  same  model  was  used  for  Winter,  Summer  and  Tropical  conditions). 
This  is  certainly  not  correct  since  the  aerosol  size  distributions  and 
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Figure  19.  Calculated  SWIR  Spectral  Atmospheric  Transmission  for  Various  Atmospheric  Models; 

8.5  km  Visibility,  0.75  km  Altitude,  5 km  Range 


Figure  21.  Calculated  SWIR  Spectral  Atmospheric  Transmission  for  Various  Water  Continuum  Models 
Tropical  Model  Atmosphere,  8.5  km  Visibility,  0.75  km  Altitude, 

5 km  Range 
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Figure  22.  Calculated  LWIR  Spectral  Atmospheric  Transmission  for  Various  Water  Continuum  Models; 
Tropical  Model  Atmosphere,  8.5  km  Visibility,  0.75  km  Altitude, 

5 km  Range 


Figure  24.  Calculated  LWIR  Spectral  Atmospheric  Transmission  as  a Function  of  Visibility; 
Midlatitude\Summer  Model  Atmosphere,  5 km  Range,  0.75  km  Altitude 
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optical  properties  are  known  to  vary  with  these  conditions.  However, 
since  aerosol  attenuation  will  be  shown  to  have  little  effect  on  relative 
performance  for  the  Lowtran  model  used  here,  this  simplification  will 
not  introduce  any  significant  additional  error.  If  a more  sophisticated 
aerosol  model  were  to  be  used  however,  this  inconsistency  should  be 
corrected. 


3.  Background  Radiance 


The  denominator  of  Equation  4-2,  describing  Barhydt's  performance 
model  involves  the  radiance  of  the  background.  Background  radiance, 
for  a thermal  imager  viewing  ground  targets  through  a long  atmospheric 
path,  is  comprised  of  radiance  from  the  terrain  (L  ),  and  from  the  at- 

o 

mospheric  path  (L  ).  The  former  component  is  attenuated  by  the  at- 
mosphere,  while  the  latter  component  (as  it  is  generally  defined  and 
calculated*)  is  not: 


(4-10) 


Previous  analyses  of  this  type  (see  Section  III),  have  neglected  the  path 
radiance  term,  and  apparently  also  neglected  the  atmospheric  attenua- 
tion of  the  terrain  radiance: 

Although  these  two  approximations  are  compensating,  the  end  result 
certainly  is  not  exact.  The  effect  of  this  is  discussed  in  Section  V. 


i 


a.  Terrain  Radiances 

Previous  analyses  of  this  type  have  assumed  a simple  300  K black- 
body  function  for  the  terrain  (and  with  the  above  approximations,  the 

* Path  radiance  is  calculated  by  considering  simultaneous  absorption 
and  emission  of  each  segment  of  the  path.  Thus,  the  attenuation  is 
accounted  for  in  the  calculation. 
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background  radiance).  Although  this  is  certainly  a good  approximation 
for  first  order  analyses  (especially  at  long  wavelengths),  it  can  deviate 
significantly  from  reality.  In  particular,  there  can  be  considerable 
solar  reflection  in  the  short  wavelength  region  during  the  day.  In  the 
present  study,  actual  measurements  of  terrain  spectral  radiant  inten- 
sities by  Smith  and  Blay  [59]  have  been  employed.  These  measure- 
ments were  taken  from  relatively  short  ranges  (~  80m),  so  atmospheric 
attenuation  effects  were  minimal.  Both  day  and  night  measurements 
were  made  with  two  high  resolution  interferometer-spectrometers  for 
the  two  spectral  regions  (see  Figure  25). 
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b.  Daytime  Measurements 


Figure  25.  Measured  Spectral  Radiant  Intensity  of 
Typical  Terrain  Backgrounds  by  Smith  and  Blay  [59], 
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These  measurements  were  converted  to  radiances  in  quantum 
units  by  the  following  procedure.  First,  the  measurements  were 
normalized  to  their  effective  radiating  areas  by  dividing  by  the  field 
of  view  - range  squared  product  to  give  radiance: 


/ <-2> 

These  results  were  then  roughly  fit  with  blackbody  curves  of  tempera- 
ture T and  the  effective  spectral  emittance  of  the  backgrouhd  de- 
eff’ 

rived: 


c(X)  = 


/ 


Teff} 

t>X 


BB 


(4-13) 


Using  the  effective  temperature  of  the  blackbody  function  (Teff),  and 
the  derived  emittance  (c),  the  spectral  radiance  in  quantum  units  was 

evaluated: 


(^r)g  = t(x> 


6L*(X,Teff) 


t>X 


BB 


(4-14) 


where  the  standard  Planck  function  for  spectral  radiance  in  quantum 
units  was  used: 


&L*(X,Tef{)- 

- c*x^ 

Lj  JiA . il 

6X 

BB  L1 

HA J 

Results  of  this  calculation  are  shown  in  Figure  26  for  day  and 
night  conditions  in  the  two  spectral  regions.  In  the  long  wavelength 
region,  the  two  curves  are  very  similar  as  might  be  expected,  since 
the  air  temperature  was  only  4 K cooler  during  the  night  measure- 
ment than  during  the  day  measurement.  In  the  short  wavelength  re- 
gion, the  radiances  during  the  day  are  typically  almost  a factor  of 
two  above  those  during  the  night.  Unfortunately  these  measurements 
are  only  for  one  locale  (and  one  type  of  terrain)  and  one  season. 
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Aberdeen  Proving  Grounds,  Maryland  in  September.  Thus,  arbitrarily 
using  these  values  with  atmospheric  models  for  Summer,  Winter  and 
Tropical  conditions  is  somewhat  of  an  inconsistency.  However,  it  will 
be  shown  (see  Section  V)that  the  relative  performance  of  thermal  im- 
aging systems  is  not  strongly  influenced  by  the  level  of  the  terrain 
radiance,  so  it  is  felt  that  this  inconsistency  is  not  significant. 

In  addition  to  these  measurement -derived  terrain  radiances,  theo- 
retical radiances  were  also  employed  in  some  of  the  performance  cal- 
culations. These  theoretical  radiances  were  calculated  directly  from 
the  Planck  function  (Equation  4-15)  for  various  temperatures  (see 
Figure  27). 

b.  Path  Radiances 

The  effective  radiance  for  the  atmosphere  in  the  viewing  path, 
arises  from  the  simultaneous  absorption  and  emission  of  the  consecu- 
tive portions  of  the  path  as  well  as  scattering  of  radiation  (from  sources 
outside  the  line  of  sight)  by  the  atmosphere  into  the  line  of  sight.  Cal- 
culation of  this  radiance  is  very  similar  in  many  respects  to  the  cal- 
culation of  atmospheric  transmission,  in  that  the  same  band  models 
and  aerosol  scattering  model  can  be  used.  Indeed,  an  SAI  atmospheric 
transmission  code  "Radst"  [60],  calculates  both  transmission  and  path 
radiance  simultaneously  (unfortunately  Lowtran  does  not  currently  in- 
clude this  feature  although  there  apparently  are  plans  to  incorporate  it 
in  a future  modification). 

The  computer  code  Radst  was  employed  in  this  study  to  evaluate 
the  path  radiance  contribution  to  the  background.  Modifications  to 
this  code  were  made,  so  that  the  atmospheric  and  aerosol  models  were 
consistent  with  those  of  Lowtran.  One  exception  to  this  was  the  short 
wavelength  water  continuum,  which  was  incorporated  in  the  Lowtran 
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results,  but  unfortunately  not  in  those  of  Radst.  It  was  felt  that  this 
modification  was  beyond  the  scope  of  the  present  effort,  and  would 
have  minor  impact  on  the  performance  calculations.  The  Radst  re- 
sults were  converted  to  quantum  units  (i.  e. , from  watts  to  photons/ 
sec)  in  the  subsequent  special  purpose  code  "Path.  " 

Sample  results  from  these  codes  are  shown  in  Figures  28-37 
for  various  conditions.  Figures  28-33  present  the  calculated  spec- 
tral path  radiance  as  a function  of  range  for  both  spectral  regions  and 
all  seasons.  These  results  are  similar  to  the  atmospheric  transmis- 
sion calculated  by  Lowtran  (see  Figures  13-18),  but  with  roughly  an 
inverse  character.  That  is,  where  the  transmission  is  large  the  path 
radiance  is  small,  and  vice  versa.  Also  the  path  radiance  increases 
with  increasing  range,  decreasing  visibility  and  increasing  humidity 
while  the  transmission  has  the  opposite  dependency  on  these  param- 
eters. It  may  also  be  noted  that  the  calculated  path  radiance  has  more 
spectral  structure  than  the  calculated  transmission.  This  is  an  arti- 
fact of  the  respective  codes  used  to  generate  the  results;  Lowtran  has 
lower  spectral  resolution  than  Radst.  The  upper  limits  to  the  path  ra- 
diance at  long  ranges  (e.  g.  , between  4.2  and  4.6  pm)  and  in  the  long 
wavelength  window)  corresponds  to  gray  body  curves  at  some  effective 
atmospheric  temperature.  Note  that  the  Tropical,  long  wavelength 
region  (see  Figure  33),  attains  this  condition  more  completely  than 
th  short  wavelength  region,  due  to  its  larger  opacity.  The  Winter 
long  wavelength  region  (see  Figure  29)  does  not  show  this  "blackening" 
effect  because  of  the  reduced  water  content.  This  Winter  radiances 
are  reduced  by  more  than  a factor  of  two  from  those  of  the  Summer, 
due  to  reduced  temperature  (and  for  the  long  wavelength  region,  re- 
duced water). 

Aerosol  content  (i.  e. , visibility)  has  very  little  effect  on  path  ra- 
diance (see  Figures  34  and  35).  This  can  be  explained  by  realizing  that 
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Figure  28.  Calculated  SWIR  Spectral  Atmospheric  Path  Radiance 

as  a Function  of  Range; 

Midlatitude  Winter  Model  Atmosphere 
8.5  km  Visibility,  0.75  km  Altitude 
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Figure  29.  Calculated  LWIR  Spectral  Atmospheric  Path  Radiance 

as  a Function  of  Range: 

Midlatitude  Winter  Model  Atmosphere 
8.5  km  Visibility,  0.75  km  Altitude 
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Figure  31.  Calculated  LWIR  Spectral  Atmospheric  Path  Radiance 

as  a Function  of  Range; 

Midlatitude  Summer  Model  Atmosphere 
8.5  km  Visibility,  0.75  km  Altitude 
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Figure  32.  Calculated  SWIR  Spectral  Atmospheric  Path  Radiance 

as  a Function  of  Range; 

Tropical  Model  Atmosphere,  8.5  km  Visibility, 

0.75  km  Altitude 
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Figure  33.  Calculated  LWIR  Spectral  Atmospheric  Path  Radiance 

as  a Function  of  Range; 

Tropical  Model  Atmosphere,  8.5  km  Visibility, 

0.75  km  Altitude 
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Figure  34.  Calculated  Spectral  Atmospheric  Path  Radiance 
as  a Function  of  Visibility; 

Midlatitude  Summer  Model  Atmosphere 
5 km  Range,  0.75  km  Altitude 
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aerosol  extinction  is  comprised  of  both  scattering  and  absorption. 

For  path  radiances  only  the  absorption  component,  which  is  by  far 
the  smaller  of  the  two  components,  is  considered  in  the  Radst  code. 
The  larger  scattering  component  would  theoretically  affect  the  path 
radiance  by  scattering  radiation  from  sources  which  would  normally 
be  outside  the  line  of  sight,  into  the  line  of  sight;  however  this  mechan- 
ism is  not  included  in  Radst.  Neglecting  this  phenomena  is  not  ex- 
pected to  have  much  effect  on  the  final  performance  results,  except 
perhaps  at  short  wavelengths  during  the  daytime  when  the  scattered 
sunlight  could  be  important. 

4.  Target  Signatures 

The  numerator  of  Equation  4-1,  describing  Barhydt's  perform- 
ance model  involves  the  signature  of  the  target.  The  temperature 
derivative  of  the  spectral  radiant  intensity  of  the  target  in  quantum 
units  is  the  actual  quantity  required.  Previous  analyses  of  this  na- 
ture (see  Section  in),  simply  used  the  temperature  derivative  of  the 
Planck  function  generally  corresponding  to  a 300  K blackbody.  The 
present  analysis  employed  effective  temperatures  and  projected  areas 
of  a typical  tank  (see  Table  3),  as  measured  by  Moulton  [61]. 

To  evaluate  the  function  required,  Equation  4-15  was  differen- 
tiated with  respect  to  temperature: 


t>  6L  ^»TefP  C1C2  Z u r,-  _ 

6T  m = IT  2 — ’ where  z = exP 

&TI  6X  TV  ( z^_x) 
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Then  the  resulting  expression  was  multiplied  by  the  various  projected 
areas,  and  summed  over  the  target: 


6 Tt>A 


(A  ). 

v p'l 


b TbA 


Table  3. 

Measured  Effective  Temperatures  and 
Projected  Areas  of  an  Operating  Tank  [61] 

Aspect 

Component 

Projected  Area 

Effective 
Temperature 
Above  Ambient 

<m2) 

(K) 

Side 

Turret 

3.4 

0.66 

Tread 

8.4 

5.  40 

Front 

Turret 

3.9 

0.  92 

Hull 

1.9 

1.90 

Tread 

1.8 

5.70 

Rear 

Turret 

2.6 

0.  91 

Hull 

1.8 

2.  20 

Tread 

1.8 

5.  50 

Engine 

1.4 

49.0 

This  radiant  intensity  function  was  calculated  for  three  different  as- 
pects of  the  tank  in  the  two  spectral  regions  of  interest.  It  was  tacitly 
assumed  that  the  emittance  of  the  target  is  the  same  in  the  two  at- 
mospheric windows. 

Results  of  these  calculations  for  front,  side  and  rear  aspects 
are  shown  in  Figure  36.  In  addition  to  these  calculations  using  the 
measured  temperatures,  calculations  were  performed  for  a range 
of  assumed  temperatures  for  the  same  projected  areas  (actually  the 
projected  area  is  not  important  since  it  can  be  canceled  out  of  the 
relative  performance  expression  (i.  e. , it  is  independent  of  wave- 
length). These  results  are  shown  in  Figure  37. 
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Figure  36.  Target  Spectral  Radiation  Function  (phot/sec/sr/^ni/K) 
for  Tank  at  Various  Aspects  Based  on  Measurements  by  Moulton  [61  ]- 
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Figure  37.  Calculated  Target  Spectral  Radiation  Function 
(phot/sec/sr/^m/K)  for  Various  Target  Temperatures 
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SECTION  V 

CALCULATED  RELATIVE  PERFORMANCES 

Using  the  procedures  and  models  described  in  Section  IV,  the  per- 
formance of  a thermal  imaging  system  operating  in  the  short  wave- 
length atmospheric  window,  relative  to  one  operating  in  the  long  wave- 
length window  was  evaluated  for  a number  of  conditions.  All  relative 
performance  calculations  use  a relative  detector  parameter  (see  Sec- 
tion IV,  #1),  of: 

4 = d q SWIR  = 2.83  (5-1) 

KV\  >LW1R 


It  should  be  emphasized  that  this  value  corresponds  to  advanced 
systems  postulated  by  Hughes  [19].  More  conservative  estimates 
(which  some  might  prefer)  would  simply  shift  the  relative  perform- 
ance curves  shown  later  in  this  section,  down  by  the  same  factor. 

Figure  38  shows  the  relative  performance  as  a function  of  range 
through  the  three  model  atmospheres  against  a nominal  target  (tank, 
rear  aspect)  and  a nominal  terrain  (night  time).  It  is  seen  that  the 
performance  of  the  short  wavelength  system  rapidly  improves, 
relative  to  the  long  wavelength  system,  with  increasing  range  for 
moderate  to  high  humidities  (i.  e. , Summer  and  Tropical  conditions, 
respectively).  From  zero  to  twenty  kilometers  range,  the  relative 
performance  increases  about  two  orders  of  magnitude  in  the  Tropical 
atmosphere,  and  about  one  order  of  magnitude  in  the  Midlatitude  Sum- 
mer atmosphere.  On  the  other  hand,  the  relative  performance  de- 
creases slightly  with  range  for  the  Winter  model  atmosphere  (less 
than  50%  from  zero  to  twenty  kilometers  range).  Cooler  targets  (or 
other  than  rear  aspect)  would  reduce  the  relative  performance  as  will 


be  shown  later. 


< 


Figure  38.  Calculated  Performance  of  a 3.4  - 4.8  FLIR 
Relative  to  an  8.1  - 12.2  FLIR  as  a Function  of 
Atmospheric  and  Aerosol  Models; 

0.75  km  Altitude,  Nominal  Target  and  Terrain 
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The  linearity  of  these  relative  performance  carves  (especially  at 
high  humidities),  indicates  an  exponential  relationship  with  range. 
This  is  not  unexpected,  since  most  of  the  transmission  terms  vary  ex- 
ponentially with  range.  This  exponential  behavior  applies  to  all  con- 
tinuum and  aerosol  attenuation,  and  "weak"  molecular  lines.  Molecu- 
lar lines  in  the  strong  line  region  vary  according  to  the  exponential  of 
the  square  root  of  range.  It  may  also  be  noted  that  the  curves  tend  to 
converge  at  small  ranges.  This  is  also  expected,  since  the  effects  of 
atmosphere  diminish  with  range.  The  value  at  zero  range  is  simply 
a function  of  the  target  and  terrain  temperatures  and  emissivities, 
and  the  relative  detector  parameter,  £. 

The  effect  of  aerosols  (i.  e. , visibility)  on  the  relative  perform- 
ance is  quite  small  (see  Figure  38).  This  result  was  surprising  at 
first,  but  is  felt  to  be  a consequence  of  the  approximations  used  in 
the  Lowtran  aerosol  model  (see  Section  IV,  #2  ).  That  is,  since  the 
spectral  shape  of  the  aerosol  transmission  is  invariant  with  visibility, 
the  ratio  of  the  average  short -to-long  wavelength  transmission  does 
not  change  much  with  visibility  (see  Figure  12): 

^SWIR^LWIR  = 0,96  at  23  km  visibility 

(5-2) 

Tswn/TL\viR  = 0,78  at  5 km  visibmty 

Since  the  performance  varies  approximately  with  the  square  root  of 
the  average  aerosol  transmission,  one  would  expect  a change  in  rela- 
tive performance  of  about  V 0.96/0.78  2=  10%  with  a change  in  visi- 
bility from  5 to  23  km.  This  is  on  the  order  of  what  is  observed  in 
Figure  38  at  low  relative  humidities  (where  the  square  root  relation- 
ship is  best).  It  should  be  emphasized  however,  that  this  result  is  felt 
to  be  fictitious,  and  the  result  of  an  unrealistic  aerosol  model  in  Low- 
tran. 


The  level  of  the  water  vapor  continuum  absorption  in  the  two  at- 
mospheric window  regions  is  still  uncertain  (see  Section  IV,  #2X>).  A1 
though  the  effect  of  this  uncertainty  appears  quite  small  in  atmospheric 
transmission,  (see  Figures  21  and  22)  its  influence  on  the  relative  per- 
formance is  rather  large  (see  Figure  39).  Combinations  of  various  con 
tinua  were  selected  to  give  the  maximum  variations  (i.  e.  , the  smallest 
in  one  band  with  the  largest  in  the  other  band,  and  vice  versa).  By 
this  pairing,  with  the  Tropical  atmosphere  there  is  more  than  a factor 
of  five  difference  between  the  maximum  and  minimum  relative  per- 
formance at  twenty  kilometers  range.  It  is  unlikely  that  the  calcula- 
tions would  be  off  by  this  much,  but  these  results  are  presented  to 
show  a worst  case.  Note  that  the  nominal  model  used  in  this  study 
does  not  correspond  to  any  of  those  indicated  in  Figure  39,  but  is 
near  Curve  A.  Its  closeness  to  Curve  A can  be  explained  by  the  fact 
that  the  long  wavelength  continuum  is  more  important  than  the  short, 
and  that  the  Lowtran  in  model  (which  is  used  as  nominal)  is  so  highly 
absorbing  (see  Figure  22). 

Platform  altitude  generally  influences  the  relative  performance 
more  than  visibility  (see  Figure  40).  At  higher  humidities  (i.  e. , 
Tropical  and  Summer  models),  the  short  wavelength  system  tends  to 
become  favored  at  lower  altitudes  because  of  the  increased  water 
vapor  in  the  path.  However  for  summer  conditions,  the  improvement 
is  less  than  a factor  of  two  at  twenty  kilometers  range  for  a decrease 
in  altitude  from  1 to  0.5  km.  An  even  smaller  change  is  indicated 
for  winter  conditions,  where  the  effect  is  insignificant.  This  is  due 
to  the  reduced  water  content  of  the  atmosphere  during  the  winter.  It 
is  interesting  (but  still  not  significant),  that  the  relative  performance 
improves  slightly  with  altitude  during  the  winter. 

Figure  41  (solid  curves)  shows  the  effect  of  target  aspect  on  rela 
tive  performance  for  nominal  atmospheric,  viewing  and  terrain  con- 
ditions. The  aspect  effect  is  quite  small  for  this  target  (less  than  a 
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Figure  40.  Calculated  Performance  of  a 3.4  - 4.8  (xm  FLIR 
Relative  to  an  8.1  - 12.2  Mm  FLIR  as  a Function  of 
Platform  Altitude;  8.5  km  Visibility, 

Nominal  Target  and  Terrain 
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factor  of  two),  and  approximately  independent  of  range  out  to  twenty 
kilometers.  Arbitrarily  varying  the  target  temperature  over  large 
extremes  (see  Figure  41,  dashed  curves),  does  result  in  significant 
performance  variations;  the  short  wavelength  system  favoring  high 
temperature  targets  and  the  long  favoring  low  temperatures,  as  ex- 
pected. The  zero  range  performance  varies  with  target  temperature 
according  to  theory  (see  Appendix  C). 

The  effect  of  terrain  temperature  is  much  smaller  than  that  of 
target  temperature  (see  Figure  42).  This  is  due  to  the  square  root 
dependence  on  terrain  radiance,  and  the  combination  of  this  term  with 
path  radiance  (see  Eqn.  4-10).  Note  also  the  reduction  of  the  influence 
of  the  terrain  with  increasing  range.  This  is  because  the  path  radiance 
term  grows  with  range,  and  starts  to  dominate  the  terrain  term. 
Measured  terrain  radiances  indicate  larger  effects  from  day  to  night 
than  might  be  expected  from  their  effective  temperatures,  but  are 
still  relatively  insignificant. 

As  mentioned  earlier  (see  Section  4.3),  most  other  analyses  ap- 
proximated the  term  t L*  + L*,  by  L*  alone.  This  approximation  is 

— & P S 

seen  to  favor  slightly  the  short  wavelength  system  for  humid  condi- 
tions and  the  long  wavelength  systems  for  dry  conditions  (see  Figure 
43).  The  effect  is  quite  small  (i.  e. , 10-20%),  and  appears  independ- 
ent of  range.  On  the  other  hand,  neglecting  L*,  but  including  the  r 

P & 

factor,  results  in  an  approximation  which  for  humid  conditions,  tends 
to  favor  the  long  wavelength  system  more  and  more  as  range  increases. 
For  dry  conditions,  it  favors  the  short  wavelength  system  slightly,  and 
again  apparently  independent  of  range. 

A cursory  examination  was  also  made  of  the  influence  of  the  limits 
of  the  spectral  bandpass  on  the  performance  of  the  systems.  Four 
smaller  bandpasses  from  within  the  nominal  short  wavelength  window 
were  selected  for  investigation:  3.4  - 4.0  pm,  3.45  - 3.5  pm,  4.5  - 
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Figure  42.  Calculated  Performance  of  a 3.4  - 4.  8 ym  FLIR 
Relative  to  an  8. 1 - 12.  2 4m  FLIR  as  a Function  of 
the  Terrain;  Midlatitude  Summer,  0.75  km  Altitude, 

8.5  km  Visibility,  Nominal  Target 
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Figure  43.  Calculated  Performance  of  a 3.4  - 4.  8 ym  FLI'R 
Relative  to  an  8. 1 - 12.  2 urn  FLIR  for  Various 


4.65  pm  and  4.6  - 4.8  |im.  The  basis  of  these  regions  was  the  calcu- 
lated atmospheric  transmission  characteristics  (see  Figure  19).  Each 
narrow  spectral  band  was  selected  to  be  highly  transmitting  and  least 
affected  by  water  vapor  in  the  atmosphere. 

Figure  44  indicates  the  relative  performance  for  each  of  these 
bands  (solid  curves),  and  for  the  nominal  band  (dashed  curve).  It  is 
seen  that  for  this  nominal  viewing  condition  at  ranges  less  than  20  km, 
all  the  narrower  bands  are  inferior  to  the  nominal  band,  except  band  A 
beyond  about  14  km.  By  narrowing  the  bands,  the  integrated  energy 
on  the  detector  is  reduced.  Thus,  generally*  the  narrower  the  band 
the  lower  the  performance  at  short  ranges.  However  at  longer  ranges, 
the  better  transmission  characteristics  of  these  bands  will  begin  to 
dominate,  and  eventually  at  some  large  range  their  curves  are  all 
expected  to  cross  above  that  of  the  nominal  spectral  band. 

The  situation  in  the  long  wavelength  window  is  not  so  clear-cut, 
in  that  there  are  no  spectral  bands  which  are  very  clear  of  water  vapor 
attenuation  (see  Figure  20).  Thus  the  narrower  bands  for  this  region 
were  selected  arbitrarily  for  a parametric  type  analysis  (see  Figure 
45).  Here  since  the  long  wavelength  system  is  being  optimized,  the 
lower  curves  indicate  better  performance.  It  is  seen  that  by  shorten- 
ing the  long  wavelength  cutoff  (e.  g. , Curves  G and  H),  some  improve- 
ment over  the  nominal  system  may  be  realized  at  longer  ranges.  Even 
Curve  E (and  possibly  D)  appear  to  be  converging  toward  the  nominal 
curve.  However  lengthening  the  short  wavelength  cutoff  (e.  g. , Curves 
F,  B and  A)  appears  to  degrade  the  system  even  more  at  longer  ranges 
(i.  e. , the  curves  seem  to  diverge  with  increasing  range). 

The  present  results  (from  Figure  38)  are  compared  with  the  re- 
sults of  previous  analyses  in  Figures  46  and  47.  It  should  be  remem- 
bered that  the  present  results  are  for  slant  paths  from  0.75  km  altitude, 

* The  exception  (i.  e.  , Curve  B crossing  Curve  A)  is  due  to 
the  shape  of  the  Planckian  function. 
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Figure  44.  Calculated  Performance  of  a 3.4  - 4.8  pm  FLIR 
Relative  to  Several  Long  Wavelength  FLIRs; 
Midlatitude  Summer  Atmosphere,  8.5  km  Visibility, 

0.75  km  Altitude,  Nominal  Target  and  Terrain 
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Figure  45.  Calculated  Performance  of  Several  Short  Wavelength 
FLIRs  Relative  to  an  8. 1 -12.2/m  FLIR; 

Midlatitude  Summer  Atmosphere,  8.5  km  Visibility, 

0.  75  km  Altitude,  Nominal  Target  and  Terrain 
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while  the  other  analyses  are  for  sea  level  conditions.  Extrapolating 
the  results  shown  in  Figure  39,  indicates  that  sea  level  conditions 
favor  the  short  wavelength  systems.  Mundie’s  calculations  (triangles 
in  Figure  46),  are  for  contrast  and  detectability  criterion  (upper  and 
lower  triangles,  respectively),  as  expressed  in  Equation  3-2.  These 
values  are  taken  from  Mundie's  Figure  12  [2],  corresponding  to  the 
nominal  target  temperature  used  here,  and  scaled  up  by  the  factor  i. 
The  range  was  not  explicitly  stated  in  Reference  2,  but  was  taken  to 
be  one  mile  (1.6  km),  since  that  seemed  to  be  his  nominal  range.  It 
is  interesting  that  these  two  criteria  give  such  divergent  results  (al- 
most an  order  of  magnitude  different). 

Barhydt's  results  (see  Table  2)  are  also  shown  in  Figure  46  (dia- 
monds and  squares  corresponding  to  aperture-resolution  products  of 
0.8  and  1.25  in-mrad,  respectively).  It  is  difficult  to  compare  these 
results  directly,  since  they  are  for  different  atmospheric  conditions 
than  the  present  calculations.  These  values  seem  to  indicate  consider- 
ably poorer  relative  performance  than  most  of  the  other  results,  par- 
ticularly at  higher  humidities.  This  is  probably  due  to  a combination 
of  Barhydt's  lower  target  temperature,  and  omission  of  the  long  wave- 
length water  continuum.  * 

Schnitzler's  results  are  indicated  in  Figure  46  by  the  dashed  lines 
for  the  standard  Tropical  Atmosphere  at  two  visibilities.  These  curves 
were  interpreted  from  Schnitzler’s  Figure  10  for  a 3.1m  target,  by 
taking  the  inverse  ratio  of  AT_n  (see  Equation  3-4).  Reasonably  good 
agreement  is  evident  with  the  present  results  at  intermediate  ranges. 
The  slopes  of  Schnitzler's  curves  are  considerably  larger  than  those 
of  the  present  study.  This  is  attributed  to  several  effects.  First,  his 

* Yates  and  Taylor’s  measurements  (employed  by  Barhydt  and 

Mundie)  were  reduced  in  such  a way  to  omit  the  long  wave- 
length water  continuum  [17,  27]. 
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calculations  are  for  the  higher  effective  humidity  at  sea  level  as  com- 
pared to  the  slant  path  considered  in  the  present  study.  Second,  the 
Lowtran  II  model  he  used  neglects  the  4 pm  water  continuum,  while 
the  present  results  include  it.  Finally,  his  results  are  for  a slightly 
different  short  wavelength  bandpass  (3.33  - 4.17  pm  instead  of  3.4  - 
4.8  pm),  which  tends  to  give  a larger  slope  (see  Figure  44).  The  zero 
range  results  are  somewhat  lower  than  any  of  the  other  results,  prob- 
ably due  to  different  assumed  systems  parameters,  and  a lower  target 
temperature. 

Milton’s  results  (see  Figure  47)  also  indicate  a slightly  steeper 
slope  (at  the  high  humidity  condition)  than  the  present  results,  but  not 
as  steep  as  Schnitzler's.  The  reasons  are  probably  the  same  as  the 
first  two  given  above  for  Schnitzler.  Otherwise  the  agreement  between 
Milton's  results  and  the  present  results  is  very  good  at  the  high  humid- 
ity conditions.  At  low  humidities  Milton's  results  are  somewhat  below 
the  present  results,  but  indicate  the  same  general  trend.  The  effect  of 
aerosols  (i.  e.  , visibility)  is  much  stronger  in  both  Schnitzler's  and 
Milton's  results,  than  in  the  present  results.  The  reason  for  this  is 
unclear,  but  may  be  due  to  differences  in  the  way  atmospheric  trans- 
mission is  applied  in  the  models  (see  Equations  4-10  and  4-11),  or  dif- 
ferences between  the  Lowtran  II  and  III  aerosol  models. 
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Figure  47.  Comparison  of  Calculated  Relative 
Performance  with  that  of  Milton 
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SECTION  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 


As  shown  in  the  previous  section,  the  short  wavelength  infrared 
thermal  imaging  systems  appear  to  offer  some  advantages  over  long 
wavelength  systems  for  some  situations.  In  particular,  high  humidity, 
long  path  length,  low  altitude  operations  favor  the  short  wavelength 
system  for  the  nominal  target  considered  here  (tank,  rear  aspect). 
Target  temperature  can  have  a significant  effect  on  relative  per- 
formance, while  terrain  temperatures  are  indicated  to  have  minimal 
effect.  Although  aerosols  are  not  predicted  to  have  much  influence 
on  relative  performance,  this  is  thought  to  be  a ficticious  result  of 
the  simplified  aerosol  model  used  in  Lowtran.  A brief  investigation 
of  spectral  bandpass  optimization  showed  there  could  be  some  im- 
provement in  both  bands,  but  nothing  very  significant  for  ranges 
less  than  twenty  kilometers. 

Comparisons  with  previous  analyses  indicate  general  agreement, 
considering  the  differences  in  conditions  assumed.  The  main  discrep- 
ancy, the  slightly  larger  rate  of  increase  in  relative  performance  with 
range,  is  attributed  to  the  higher  humidity  sea  level  conditions.  Once 
again  it  should  be  emphasized  that  analyses  of  this  type  do  not  address 
the  absolute  performance  of  these  systems,  only  relative  performance. 
Thus,  the  results  presented  here  should  not  be  taken  to  indicate  that 
either  system  would  necessarily  be  capable  of  operating  for  all  the  con- 
ditions considered.  On  the  other  hand,  even  though  the  short  wavelength 
system  may  be  slightly  inferior  to  the  long  wavelength  system  at  some 
conditions,  it  may  still  have  adequate  performance  there.  This  is  quite 
likely  the  case  for  the  Winter  condition,  since  the  performance  of  both 
systems  should  be  excellent  in  this  low  humidity  environment. 


Based  on  the  above  conclusions,  it  is  recommended  that  future 
studies  of  this  type  employ  a more  sophisticated  and  better  verified 
atmospheric  transmission  and  path  radiance  model,  especially  with 
regard  to  aerosol  and  water  continuum  modeling.  In  addition,  it  may 
be  worthwhile  to  investigate,  in  more  depth,  the  selection  of  optimum 
spectral  operating  regions,  possibly  considering  several  rather  narrow 
bandpass  filters.  A statistical  analysis  of  the  effects  of  meteorological 
conditions  for  various  likely  operational  theaters  would  also  seem  to  be 
of  value  for  the  optimum  design  and  implementation  of  thermal  imaging 
systems. 
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APPENDIX  A 

AVAILABILITY  OF  INTERMEDIATE  RESULTS 


h 


5 


i. 


The  analysis  presented  in  this  report  can  be  considered  in  two 
parts:  1)  the  evaluation  of  the  inputs  to  the  performance  model:  at- 
mospheric transmission  and  path  radiance,  and  target  and  terrain 
radiance;  and  2)  the  evaluation  of  the  performance  using  the  above 
inputs  in  Barhydt’s  model.  In  the  future  it  may  be  desirable  to  utilize 
a more  sophisticated  or  specific  (i.  e. , for  a particular  system)  per- 
formance model  with  the  results  of  Part  1,  or  simply  use  segments 
of  it  for  some  other  purpose.  For  these  reasons  the  results  of  Part  1 
have  been  saved  on  magnetic  tape  for  possible  future  use. 

Atmospheric  spectral  transmission  calculated  by  Lowtran  III  has 
been  stored  (in  the  order  shown  in  Table  Al)  in  files  TRANS1  and 


Table  Al. 

Calculated  Atmospheric  Spectral  Transmission 

and  Path  Radiance  Cases 
Platform 

Slant 

Visibility 

Altitude 

Range 

(km) 

(km) 

(km) 

5.0 

0.  75 

1,5, 10,  20 

8.5 

0.  75 

1,5, 10,20 

23.0 

0.  75 

1,  5, 10,  20 

8.5 

0.  50 

1,5,10,  20 

8.5 

1.  00 

1,  5, 10,  20 
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TRANS2  (Transmission)  for  the  Midlatitude  Summer  model  for  the 
short  and  long  wavelength  atmospheric  windows,  respectively.  Simi- 
larly, files  TRANS3  and  TRANS4  contain  the  calculations  for  the  Mid- 
latitude Winter  model,  and  TRANS5  and  TRANS6  for  Tropical.  The 
short  wavelength  transmissions  are  digitized  every  0.01  |im  from  3.4 
to  4.8  pm,  while  the  long  wavelength  ones  are  every  0.02  pm  from  8.1 
to  12.2  pm.  In  order  to  adjust  the  long  wavelength  continuum,  two 
records  are  provided  for  each  case  in  TRANS2,  TRANS4  and  TRANS6: 
the  first  is  the  total  transmission  divided  by  the  transmission  due  to 
continuum  alone,  and  the  second  is  the  total  transmission. 

Atmospheric  spectral  path  radiance  calculated  by  Radst  (and  modi- 
fied by  Path)  has  been  stored  in  the  same  order  in  files  PRAD1-PRAD6 
(Path  Radiance).  The  only  difference  is  that  there  is  only  one  record 
for  each  case  for  both  the  short  and  long  wavelength  windows.  Each 
record  for  both  transmission  and  path  radiance  is  preceded  by  an 
identification  line  which  completely  specifies  the  case.  The  number 
of  entries  in  the  record  is  also  indicated  on  this  line. 

Terrain  spectral  radiance  calculated  by  Photon  for  the  cases  out- 
lined in  Table  A2  has  been  stored  in  file  TRAD  (Terrain  Radiance). 


Table  A2. 

Calculated  Terrain  and  Target 

Radiance  Cases 

Source 

Condition 

Temperature  (K) 

Terrain 

Night  and  Day 

— 

Terrain 

Theoretical 

250,  270,  290,  310 

Target 

Side  and  Front 
Rear  Aspects 

— 

Target 

Theoretical 

250,  275,  325,  350 
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The  first  record  contains  the  wavelengths  (in  pm)  for  which  these  data 

(and  the  other  data)  have  been  digitized.  The  long  wavelength  region 

is  merged  into  the  end  of  the  short  wavelength  region  in  the  same  record. 

The  second  and  third  records  contain  the  day  and  night  terrain  radiances 
o 

(phot/sec/sr/m  /pm)  as  measured  by  Smith  and  Blay  [59].  The  fourth 
through  seventh  records  have  the  calculated  blackbody  radiances  for 
temperatures  of  250,  270,  290  and  31  OK,  respectively. 

The  target  spectral  radiation  function  > Phot/sec/sr/Mm/K^ 

as  calculated  by  Photon  for  the  cases  outlined  in  Table  A2  has  been 

digitized  in  the  file  TARSIG  (Target  Signature).  The  first  three  records 

contain  the  radiation  function  for  the  typical  tank  at  front,  side  and  rear 

aspects,  respectively,  based  on  the  measurements  of  Moulton  [61]. 

The  last  four  records  have  the  radiation  function  for  a blackbody  of 

2 

projected  area  15.7  m at  temperatures  of  250,  275,  325  and  350  K, 
respectively. 
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APPENDIX  B 

CALCULA TIONA L PROCEDURE 


The  general  calculations!  procedure  is  outlined  in  the  flow  dia- 
gram shown  in  Figure  Bl.  A brief  description  of  each  of  the  com- 
puter codes  (depicted  by  circles  and  ovals  in  Figure  Bl),  is  given 

below: 

1.  LOWTRAN  HI  - The  third  generation  AFCRL  code  to  calculate 
atmospheric  spectral  transmission.  Inputs  are  the  atmospheric 
parameters  (e.  g. , seasonal  and  latitudinal  model,  and  visibility), 
viewer  slant  path  (altitude,  range),  and  spectral  bands. 

2.  RADST  - A code  developed  by  Anding  to  calculate  the  atmospheric 
spectral  path  radiance.  Inputs  are  the  same  as  those  of  Lowtran 
ID  (see  #1,  above). 

3.  PHOTON  - Calculates  the  radiant  intensity  function  of  general  tar- 
gets, and  the  radiances  of  backgrounds  from  their  temperatures 
and  emissivities  using  the  Planck  function.  Also  produces  plots  as 
a function  oi  wavelength. 

4.  TRAN  - Re-formats  the  Lowtran  output  to  even  wavelength  intervals 
and  adds  the  3-5  pm  water  continuum  (or  alternate  8-12  pm  water  con- 
tinua).  Also  produces  plots  of  atmospheric  spectral  transmission. 

5.  PATH  - Converts  RADST  output  to  quantum  units,  and  re-formats  it 
to  even  wavelength  intervals.  Also  produces  plots  of  atmospheric 
spectral  path  radiance. 

6.  FLIPER  - Calculates  relative  FLIR  performance  for  various  spectral 
bands  (within  limits  of  bands  above).  Has  the  option  of  suppressing  path 
radiance  and/or  transmission  in  the  background  radiance  term,  or  all 
atmospheric  effects  (i.  e. , zero  range  condition). 

7.  PLOPER  - Plots  performance  as  a function  of  range  for  various  con- 
ditions. 


i 


■ 
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Figure  Bl. 

Flow  Diagram  of  Calculational  Procedure 
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APPENDIX  C 

RELATIVE  PERFORMANCE  AT  ZERO  RANGE 


At  zero  range,  the  relative  performance  of  short -to-long  wave- 
length thermal  imaging  systems  is  independent  of  the  atmosphere 
(Ta  = 1 and  L*  = 0),  so  a closed  form  expression  can  be  written: 

Po  * «Jkx/Jk*2>  <LAx/LMj)1/2 
where 


j.  » r - 

JAX  6TbX  dX 


6L(T  ,b) 

L*  — sf-dX 

AX 

These  two  modified  Planckian  functions  can  be  evaluated  using  series 
expansions 


3C1T2 


WV4  x w 

r,<.Tr  Z 


C2"  n=l 


-IH  .I1  i1  S , 2 , 

6 (X  + T + T+-J> 
n n 


C1T3  - 


1 & A _ 

L’  - f 1 V / x 2x  2 k 

L'!t7T.s'  (-r  + T + -j) 


c2"  n=i 


n n 


where 


x = °2 

x " XT 


For  €tWlR  ' €t}LWIR  and  VSWIR  = *g^LWIR’  with  Tt  = 350 


and  T = 300: 
g 


P0-1.9. 


(C-l) 

(C-2 
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